NACA RM A52F16 


SECURITY INFORMATION 



NACA 


RESEARCH MEMORANDUM 


EFFECT OF CAMBER AND TWIST ON THE STABILITY CHARACTERISTICS 

c v h 

OF MODELS HAVING A 45° SWEPT WING AS DETERMINED fey ® - 


w 


THE FREE-FALL METHOD AT TRANSONIC SPEEDS 

By Maurice D. White 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 


CLASSIFIED DOCUMENT 

This material contains information affecting the National Defense of the United States within the meaning 
of the espionage laws, Title 18, U.S.C., Secs. 793 and 794, the transmission or revelation of which in any 
manner to an unauthorized person is prohibited by law. 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

August 29, 1952 




1J 


MCA RM A52F16 


CONFIDENTIAL 


MTIOML ADVISORY COMMITTEE FOR AEROMUTICS 
RESEARCH MEMORANDUM 

EFFECT OF CAMBER AND TWIST ON THE STABILITY CHARACTERISTICS 
OF MODELS HAVING A L5° SWEPT WING AS DETERMINED BY 
THE FREE-FALL METHOD AT TRANSONIC SPEEDS 
By Maurice D. White 


SUMMARY 


Measurements were made in free fall at transonic speeds of the 
dynamic stability characteristics of three models. Two of the models 
had 45 ° swept wings of aspect ratio 6 and 45 ° swept tail surfaces, and 
differed from each other only in that the wing of one of the models was 
plane and the wing of the other was cambered and twisted. The third 
model had the same fuselage-tail arrangement as the others, but had no 
wing. Static and dynamic longitudinal stability characteristics of the 
models were determined at Mach numbers ranging from 0.84 to 1.16, wing 
Reynolds numbers ranging from 2,700,000 to about 6,000,000, and angles 
of attack from about zero lift to stalling angles at the lower Mach num- 
bers, and to a maximum of 10 to 12° at the higher Mach numbers. The 
results showed no significant difference in the stability characteristics 
of the models due to cambering and twisting the wing. There was an 
appreciable variation in static longitudinal stability of the two wing-on 
models in traversing the Mach number range, with the maximum stability 
occurring at a Mach number of about 0.97. The static directional stabil- 
ity of the three models was relatively unaffected by Mach number variation. 


INTRODUCTION 


The prediction of the dynamic stability characteristics of airplanes 
or missiles is particularly difficult in the transonic speed range where 
theoretical methods are generally inadequate, and aerodynamic derivatives 
tend to change value rapidly or erratically with changing Mach number. 
References 1, 2, and 3 present experimentally determined variations of 
dynamic stability characteristics with Mach number which illustrate the 
erratic nature of the results in this region. Since it is apparent that 
experimental results for specific configurations will have to be relied 
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upon as the basis for design for some time to come, any information that 
can be added to the existing fund of data should be of value. 

Some additional experimental data have been obtained on the dynamic 
stability characteristics of several models as part of an investigation 
that was made at Edwards Air Force Base by the Ames Laboratory of the 
NACA using the free-falling recoverable-body technique. The investiga- 
tion was made to compare the characteristics of two models having high- 
aspect-ratio 45° swept wings, one plane and the other cambered and 
twisted. In addition, the same fuselage-tail arrangement was tested 
without a wing. 

The tests covered a Mach number range from 0.84 to 1.16, and a wing 
Reynolds number range from 2,700,000 to about 6,000,000. The angles of 
attack ranged from about zero lift to stalling angles at the lower Mach 
numbers, and to a maximum of 10° to 12° at the higher Mach numbers. 


NOTATION 


a x 

a z 

a . c . 

c 

C.g. 


m 

M 


longitudinal acceleration, units of g 

vertical acceleration, units of g 

aerodynamic center position, percent TT 

wing span, feet 

local wing chord, feet 

center of gravity position, percent cf 


wing mean aerodynamic chord 



feet 


acceleration of gravity, 32.2 feet per second squared 

moment of inertia in pitch, slug-feet squared 

tail length, feet 

mass, slugs 

Mach number 
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P 

q. 

Qo 

4 


period of oscillation, seconds 

angular velocity in pitch, radians per second 

dynamic pressure (W), pounds per square foot 

|sl 

5t 


r ordinate of fuselage, inches 

S wing area, square feet 

S-t horizontal-tail area, square feet 

t time, seconds 

Tjyg time for oscillation to damp to one-half amplitude, seconds 
V velocity, feet per second 

W model weight, pounds 

x longitudinal station, inches 

y spanwise distance from plane of symmetry, feet 

a single of attack, degrees 

a^. angle of attack of tail, degrees 

a T trim angle of attack, degrees 

• da , . 

a. , radians per second 

(3 angle of sideslip, degrees 

6 horizontal tail deflection, degrees 

e angle of downwash, degrees 

p atmospheric density, slugs per cubic foot 

T}-|- tail efficiency 

Cp lift coefficient 


lift \ 
<Jo S / 
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°Lt 


lift coefficient of tail 


lift 

9o S 


^m 






-m 


a 


Um s 


Cm d 


V 


pitching-moment coefficient of complete model about model center 
'moment ' 


of gravity 


qoSc 


C mw pitching-moment coefficient of exposed wing panels 


moment 


yawing-moment coefficient of complete model about model center of 

.. /moment' 
gravity 


3Cl 


<JoSc 


lift-curve slope )> P er degree 


pitching-moment-coefficient slope [ P er degree 


SCm\ 




control-surface effectiveness parameter P er degree 

5c m 


S(qc/2V) 




d(dc/2V) 


# per radian 


9 per radian 


yawing-moment-coefficient slope [ J 9 per degree 


Sc n » 


MODEL 


Three model configurations were tested, all having the same fuse- 
lage and tail-surface dimensions. Two of the configurations had wings 
and one had no wing. Figures 1 and 2 show a three-view drawing and a 
photograph of one of the winged configurations, and table I lists the 
physical specifications of the three models. The two winged models 
differed from each other only in that the wing of the plane wing model 
was symmetrical and untwisted, while the wing of the other model had 
camber and twist. For the model having the cambered and twisted wing, 
the washout of 10° at the tip (measured streamwise) was obtained by 
twisting the wing so that the constant-percent-chord lines remained 
straight. The wing was constructed of solid aluminum alloy. 
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The fuselage was 210.5 inches in length and nad a fineness ratio 
of 12.4. The fuselage ordinates from the 8 -inch to the 139*4-inch 
station are given by the equation in figure 1. From station 139.4 the 
fuselage tapered conically to a radius of 5.2 inches at station 189 . 6 . 
From stations 189.6 to 210. 5, a tail shape approximating that given by 
the equation for values of x from I 83 .I to 204 was used. 

• 

Both of the horizontal-tail surfaces and both of the vertical-tail 
surfaces of each model were all-movable, pivoting on axes perpendicular 
to the fuselage axis. A schedule of horizontal- tail movement was preset 
so that the tails would deflect and return to trim position in rapid 
pulse-type movements at regular time intervals during the test phase of 
the drop. The vertical-tail surfaces were actuated differentially by the 
roll-position stabilization system to provide roll control. All the tail 
surfaces were constructed of solid aluminum alloy. 


INSTRUMENTATION 


NACA continuously recording flight instruments were used to record 
the various quantities measured. A listing of the quantities and of the 
instruments used to measure them follows: 

Quantity Instrument 


Angle of attack and 
angle of sideslip 


Rate of pitch and 
rate of roll 


Slave selsyns or recording oscillographs 
(depending on the installation) recording 
movements of vanes mounted on boom ahead 
of body (fig. l) 

NACA turnmeter 


Angular acceleration 
in pitch 


Angular accelerometer with recording 
oscillograph 


Vertical and longitudinal Linear accelerometers with recording oscillo- 
accelerations " “ “ 


Transverse acceleration 


graph and NACA 3-component accelerometer 
NACA 3- component accelerometer 


Horizontal- and vertical- NACA 2-component control position recorder 
tail deflections 


Mach number and dynamic 
pressure 


NACA 6 - cell manometer 
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The airspeed system was calibrated at different angles of attack using 
the SCR 584 radar installation of the NACA High-Speed Research Station 
at Edwards Air Force Base. All the flight records were synchronized by 
means of a chronometric timer. 


TESTS 


The results presented in this report were obtained during a series 
of free-fall drops of the models in which the models were trimmed at 
different angles of attack and longitudinal disturbances were produced 
by intermittent movement of the horizontal-tail surfaces. The models 
were released from a carrier airplane at an altitude of 40,000 feet and 
allowed to fall freely at about zero lift attitude until the desired 
Mach number was reached. At that time the horizontal tail was moved 
abruptly to the setting for trim and, thenceforth, was pulsed at inter- 
vals of 2.4 seconds. A typical time-history plot of the control deflec- 
tions and the resultant motions of the model for a portion of a flight 
is shown in figure 3 . The results presented herein were obtained from 
analysis of the free-oscillation characteristics following each pulse. 

At the conclusion of the test periods of each drop, the models were 
recovered by the use of a dive brake and parachute. 

The models were roll-position stabilized throughout the drop. The 
system employed stabilized the model within bank angles of approxi- 
mately 10°, and within roll rates of about 0.9 radian per second. 

The airspeed system of the models was calibrated throughout the 
test range using the NACA radar-phototheodolite method. 

For the winged models, angles of attack ranging from about zero 
lift to the stall were covered for Mach numbers up to about 0.9; at 
higher Mach numbers angles of attack up to 10° were covered for the 
plane-wing model, and up to 12° for the cambered- and twisted-wing model. 
For the wing-off model, the results were obtained for angles of attack 
ranging from about -1° to 6°. 

Generally, the angle-of -attack range covered by the oscillations 
during a drop was of the order of ±4 . For the plane-wing models the 
range of angles covered at lower Mach numbers increased to as much 
as ±7° for one of the drops. 

For the winged models a Mach number range from 0.84 to 1.10 was 
covered with Reynolds numbers ranging from 2,750,000 at the lower Mach 
number to 5 >600,000 at the higher Mach number. For the wing-off model 
the Mach numbers ranged from O .98 to l.l6 with Reynolds numbers (based 
on the wing mean aerodynamic chord) ranging from 3 > 900,000 to 6 , 250 , 000 . 
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ACCURACY 


Based on the uncertainties that are estimated to have been present 
in the various component quantities, it is believed that the accuracy of 
a single determination of any quantity is as noted below: 

a ±l/4° 


% 


±0.005 (for linear range of angles of attack) 


C m ±0.012 

a 

a.c. position ±0.167 


Cm q + Cm d 


±30 percent of actual value 


Cn p 


t 


±0.003 



±0.02 


ANALYSIS 


Following the treatment used by other investigators, the flight 
data have been analyzed under the assumption that the motions of the 
body are adequately described by a linear second-order system. A 
detailed description of the method of analysis used is given in appen- 
dix A. Reference 4 presents a fairly complete discussion of this method 
of analysis and of the assumptions involved in its use. 

For some of the data included in this report, the assumptions of 
linearity inherent in the methods of analysis have been violated. Aside 
from the fact that the data have been obtained generally under conditions 
of changing Mach number, dynamic pressure, and altitude, some of the 
oscillations encompassed ranges of angle of attack foi* which the lift 
coefficient and pitching-moment coefficient did not vary linearly with 
angle of attack. 

It is apparent that the degree to which the effects of such non- 
linearities will be evidenced in the final results will depend on the 
degree of nonlinearity present. No attempt is made in this report to 
analyze the data quantitatively in terms of the degree of nonlinearity. 
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RESULTS AND DISCUSSION 
Lift 


Figure 4 shows the variation of lift coefficient with angle of 
attack at different Mach numbers for the three model configurations with 
different control settings. A significant feature of these curves is 
the decrease in slope that is evidenced by the curves for the wing-on 
models at the higher angles of attack. This characteristic is particu- 
larly evident at the lower test Mach numbers (M < 0.96). 

In figure 5 the lift-curve slopes for small angles of attack as 
determined from the data of figure 4 are presented as a function of Mach 
number. The curves of CL a for the two wing-on models show a general 
tendency to decrease with increasing Mach number through the test range. 

At the two extremes of the Mach number range the values for the two 
wings are in agreement, and at intermediate Mach numbers the values 
of c L a for the Plane-wing model appear to be only slightly greater 
than those for the cambered- and twisted-wing model. Measurements of 
the forces on the exposed wing panels as determined from a balance 
within the fuselage show trends that essentially parallel those described 
above, which indicates that the variations are due primarily to the wings. 

For the wingless model the lift-curve slope decreased generally with 
increasing Mach number between M = 0.98 and 1.16. 


Longitudinal Trim 


The trim angles of attack for all the drops were determined as the 
mean angle of attack of the oscillations. The variations with Mach num- 
ber of the trim angles of attack as determined by this method are shown 
in figure 6 for all the drops. The results show the variations of trim 
angle of attack with Mach number to be generally small throughout the 
test range for all the models . Some increase in the variation of trim 
angle of attack with control deflection is noted at the lower Mach num- 
bers. This is due largely to the decreased static stability which 
existed at these Mach numbers, rather than to increased control effec- 
tiveness . 


Static Longitudinal Stability 


Presentation of results.- The values of the static stability param- 
eter Cn^ were determined from the half-periods of the oscillations. 
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as described in appendix A, and are shown in figure J as a function of 
Mach number for each of the configurations. In evaluating the results 
it was found that the variations of period with time following succes- 
sive control disturbances did not always form a continuous curve. The 
causes of this phenomenon are still under investigation. 

For presentation in figure 7 mean lines were drawn through all the 
data points following each control deflection and these curves were 
plotted as a function of Mach number. Each short line in figure 7 rep- 
resents the variation of C ma following a control movement. This method 
of presentation illustrates the discontinuous nature of the data which 
was described above and shows also the variation of Mach number during 
the interval following a control disturbance. A mean line drawn through 
the short curves appears to indicate some of the larger trends of the 
data, even though small changes cannot be accurately identified. 

As a result of operation of the roll stabilization system, varia- 
tions in rolling velocity occurred during the oscillations with the 
maximum roll rates reaching values as large as 0.9 radian per second. 
Calculations based on reference 5 using measured frequencies in pitch 
and yaw indicate that a steady rate of roll of 0.9 radian per second 
would affect the value of Cm a by only about 5 percent. Whether the 
effects of the oscillatory roll actually experienced would be greater 
or less remains to be established. 

The faired curves of Cj^ in figure 7 were combined with effective 
values of Cp a for the particular data to compute the aerodynamic cen- 
ter variations plotted in figure 8. 

Effect of Mach number .- For the wingless model the negative values 
of Cm a decreased progressively with increasing Mach number over the 
entire test range of Mach numbers despite the rearward movement of the 
aerodynamic center position that occurred over part of the Mach number 
range (M = 1.06 to M = 1.16). The decrease in lift-curve slope that 
occurred over this range (21 percent decrease between M = 1.06 
and M = 1.16) was apparently great enough to more than offset the effect 
of the aerodynamic center movement. 

The variations of Cm a and aerodynamic center position with Mach 
number were of the same general character for both the winged models at 
small angles of attack. The values of increased negatively with 

a corresponding rearward movement in aerodynamic center position as the 
Mach number was increased up to about M = 0.98* As the Mach number 
increased further and the aerodynamic center moved slightly forward, the 
values of Cm a decreased negatively. The differences in variation 
of C ma between the cambered- and twisted-wing and plane-wing models 
were somewhat obscured by the scatter of the data. The range of values 
of c m a covered by the two winged models, however, appeared to be 
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essentially the same and somewhat greater than that noted for the wing- 
less body over the restricted Mach number range for which comparisons 
could be made. 

Effect of angle of attack .- As may be seen from figure 6 a number 
of tests were made in which the oscillations centered around different 
angles of attack. The results in figure 7 indicate that except at Mach 
numbers below M = 0.92 (with corresponding Reynolds numbers below 
3,300,000) no large or consistent effects of angle of attack could 
be distinguished in the data. At the Mach numbers below M = 0.92 the 
data indicate a difference in the variation of Cm a with angle of attack 
for the two winged models; that is, for oscillations covering high angles 
of attack, the cambered- and twisted-wing model showed an apparent 
decrease in negative value of Cm a , while the plane-wing model showed an 
apparent increase. This difference in behavior is probably associated 
with the facts that (l) the angle-of-attack ranges designated as high and 
low were not the same for the two models, and (2) the variations of C m 
with a were nonlinear at these Mach numbers. Figure 9 shows, at two 
Mach numbers, the variation of pitching -moment coefficient with angle of 
attack for the exposed wing panel as measured by a balance within the 
fuselage. In figure 9 letters are appended to the curves for M = 0.86 
to identify the ranges of angles of attack covered in particular oscilla- 
tions, and in figure 7 the stability data associated with these angles of 
attack are denoted by the same letters; that is, the curve labelled A-B 
in figure 7(b) represents data obtained from an oscillation between the 
angles of attack indicated by the letters A and B in figure 9( a )> etc. 

It is apparent from the data of figure 9 that, whereas the oscillations 
of the cambered- and twisted-wing model covered two ranges of angles of 
attack in which the slopes of Cm w were distinctly different (E-F and 
G-H), the oscillations of the plane-wing model covered ranges of angles 
of attack over which the variation of Cm w was extremely nonlinear for 
both oscillations (A-B and C-D) . For the latter cases it would be diffi- 
cult to assign effective values of Cm^ to either the low- or the high- 
angle-of -attack oscillations, and it would therefore not be surprising 
if either an apparent increase or decrease in negative value of Cm a 
were indicated for the higher angles of attack. 

It is noteworthy that at the higher Mach numbers where no consistent 
differences due to angle of attack were evidenced by the oscillation data 
of figure 7 , the wing-panel moment-coefficient data of figure 9 also 
showed little change in slope with angle of attack. The data 
for M = 1.02 in figure 9 illustrate the relative linearity of the 
variations of Cmy- with a at higher Mach numbers. 
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Damping in Pitch 


The values of the damping-in-pitch parameter, C m ^ + C m ^, determined 
by the method outlined in appendix A, are plotted as a function of Mach 
number for the three test configurations in figure 10. For comparison 
with the experimental results, values of Cmq + Cm^ estimated as 
described in appendix B are shown in the figures. 

For the wingless body the data in figure 10(a) indicate that there 
is little variation in the value of C m + Cm^ as the Mach number is 

increased between M = O .98 and M = l.lo. The estimated values were in 
good agreement with the experimental values throughout the test range 
of Mach numbers. 

For the wing-on models figures 10(b) and 10(c) show that at Mach 
numbers above M = 1.0 the experimental values of Cm + Cm^ for the 

two models were in essential agreement, both curves showing a slight 
decrease in value with increasing Mach number. For Mach numbers less 
than M = O .96 there was a considerable difference in the level of the 
values of Cmq + Cm^ for the two models, although the variation with 
Mach number for the two models was about the same. This difference 
between the experimental values for the two models is associated with 
flight conditions for which greater nonlinearities in lift curves 
occurred for the plane-wing model than for the cambered- and twisted- 
wing model. (See fig. k.) Since the value assumed for the lift-curve 
slope directly affects the experimentally determined value of Cmq + C m< j 
(see appendix A), it is possible that the procedure employed here in 
evaluating the nonlinear lift-curve slope, that is, taking the average 
slope over the appropriate angle of attack range, gives the damping 
effect of a nonlinear lift-curve slope less weight than it deserves. 
However, even the assumption of a value of the lift-curve slope equal 
to that obtained at small angles of attack would not completely elimi- 
nate the differences between the two models. It would appear, therefore, 
that there are other differences j.n the characteristics of the models 
which are associated with the nonlinear lift curves and which affect the 
values of Cm + C m &. 

For both the wing-on models there appeared to be a localized 
decrease in the value of C m q + which occurred at a Mach number 

of O .96 for the plane-wing model and 0.98 for the cambered- and twisted- 
wing model. These Mach numbers, incidentally, correspond respectively 
to the Mach numbers where the static stability parameter Cm a reached 
a peak (fig. 7)* This local decrease in value of C m q + has been 

experienced in tests of other swept-wing models (reference 2) and, while 
of relatively small magnitude for the present models, might be of appre- 
ciable significance for airp lanes having more conventional tail volumes. 
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Comparison of the estimated values of Cm^ + Cn^ with the experi- 
mental values for the two wing-on configurations shows reasonably good 
agreement for the cambered- and twisted-wing model throughout the Mach 
number range and for the plane-wing model for Mach numbers above 
M = O. 98 . For the plane-wing model the predicted values are consider- 
ably less than the experimental values for Mach numbers less than 0.94. 
The possible reason for this latter discrepancy has already been dis- 
cussed as being at least partly due to the nonlinear lift-curve slopes. 


Control-Surface Effectiveness 


Values of the control-surface effectiveness parameter dC m /dS 
as determined by the method shown in appendix A are plotted in figure 11 
as a function of Mach number. Data are shown only for the winged models 
over a limited range of Mach numbers. Results for the wing-off model 
were not included in figure 11 because they were based on data for only 
one drop and were therefore considered too inaccurate to be compared 
with the data for the wing-on models. 

The results were consistent for the two winged models and, except 
for the decrease in value indicated at a Mach number of 0.94, the varia- 
tion with Mach number is similar to the variation in lift-curve slope 
that occurs for 45 ° swept surfaces over this Mach number range. 


Directional Stability 


In some of the tests well-defined directional oscillations were 
indicated in the sideslip-angle records. These oscillations which were 
of usable regularity only for the smaller angles of attack were analyzed 
as described in appendix A to obtain the values of the static directional 
stability parameter, Cn Q f > shown in figure 12. For comparison with these 

p 

results the values of Cn^ for the wing-off configuration from figure 7 
are also shown. The results show that the values of Cng f are muc ^ L 
smaller than the values of C m(X for the wing-off configuration. This 
would be anticipated as a result of the vertical tail being smaller in 
size than the horizontal tail. The curves of Cn^ 1 show no large varia- 
tions with Mach number for any of the configurations. 

Comparison of the values of Cnp 1 for the various configurations 
showed the cambered- and twisted-wing model to have greater directional 
stability than either of the other models for Mach numbers less than 
about 0.9. 
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CONCLUDING REMARKS 


Free-fall drop tests of three models at transonic Mach numbers 
indicated the following: 

1. There was little difference throughout the Mach number range 
in either the static or dynamic longitudinal stability characteristics 
of two 14-5° swept-wing models, one having a cambered and twisted wing 
and the other a plane wing. 

2. There was an appreciable variation in static longitudinal 
stability of the two wing-on models in traversing the Mach number range 
with the maximum stability occurring at a Mach number of about 0.97. 

3. For the test technique employed, which consisted of pulsing 
the longitudinal control at regular time intervals, the variations with 
Mach number of the periods of the short-period oscillations were not 
continuous. The causes of this behavior are not yet established. 

L. The variation of static directional stability with Mach number 
was small for the small angles of attack at which results were obtained. 
The cambered- and twisted-wing model appeared to have greater directional 
stability than the other models for Mach numbers less than about 0.9. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif. 
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APPENDIX A 

ANALYSIS OF FLIGHT DATA 
Equations of Motion 


The results in the present investigation were analyzed by assuming 
that the following equations describe the motions of the model completely: 


Iyq = MQa + Mqq + M^a + Mgb 
mV(a - q) = -L^a - Lgb 

Combining these equations, the relationship is obtained that 

q _ Cp + C]D 
5 


and where 


Co ~ 


b = 


d 2 + 

bD + k 

and 

k are d 

Mg La 

- L$Ma, 

mVIy 

Mb 

LsMd 

iy " 

mVI Y 

L a 

Mq + 

mV 

J Y 

M<x 


'lY ' 

" mVIy 


( 1 ) 

( 2 ) 

(3) 


Mb = £ 5 “ 57.3qoSc 


La = 


ST 57.3^8 


' 5? 57 - 3, 3° s = 

* - It 57 - 3, 3o s 
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M& = 


M q = 


Sc 


■m 


i ( 

\2V 


Sc 


■m 


a 'I 




o- c 

^2V 


The denominator of equation (3) defines the control-fixed oscilla- 
tory characteristics of the model, the constant b being related to the 
damping of the oscillation by the expression 


b = 


1.386 


and the constant k, to the period of the oscillation by the expression 

2rt 


P "/TU 


Evaluation of Flight Data 


The detailed procedures used in evaluating the various quantities 
discussed in this report are described in the following paragraphs. 

Lift-curve slope, CLq, .- In order to evaluate the lift-curve 
slope Cp^ time histories were made of the factors a and 




cos a 


a x sin a) 


For each of a number of values of a, Cp, was determined and plotted 
against Mach number. For selected Mach numbers, faired values of Cl 
were determined from these plots and were plotted against the corre- 
sponding angles of attack. 


Static stability parameter, Cm a . - Considering only the regions 

where the longitudinal control was held fixed, measurements were made 
of the time intervals between successive peaks of the short-period 
oscillations. These measurements were made on four different records 
(<1> 4, a z, and o-)> a nd the values were plotted as a function of time 
and were faired. Values of the dynamic pressure qQ were also plotted 
as a function of time. Finally, at sufficient points to define the 
particular curves, the values of Cma we ne computed using the relation 
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Crn = - 1 

57 - 3s ^ 0 (ir 

In the cases where the control moved despite efforts to fix it, the 
movement was found to be roughly in phase with the angle of attack and a 
correction was applied to the value of Cn^ as follows: 


AC 


SCm dS 
(38 da 


where c)C m /c3& was determined from the test data, and the values of d8/da 
from flight measurements . 


The results from tests made wirh different center-of-gravity loca- 
tions were converted to one center-of-gravity location by applying a 
correction 


AC = [(c.g.)test - (c.g.)nl 

^ 100 


CL a 


The assumptions implicit in the procedure of analyzing the data by 
the methods described above are discussed in reference 4. These include 
the assumption that the equation of motion along the longitudinal axis 
may be ignored, and that certain terms may be ignored in the equations 
with little error. The errors in Cn^ due (l) to ignoring the term (b/2) 
in the equation 


P 



and (2) to ignoring the term 


LqM q 

mVTy 


in the equation 


_ Ma 

mVIy Iy 


are at the most 3 percent and 2 percent, respectively. 


Aerodynamic center position, a.c .- The aerodynamic center position 
was calculated as a fraction of the mean aerodynamic chord by using the 
equation 


a.c. = 0.24 - 


( dCmAx A 
\bCiJ^a.y 


where dc m /ckx was obtained from the faired curves of figure 7> and the 
effective values of hCjj/^a for the appropriate angle-of-attack ranges 
were used. 
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Damping in pitch parameter, Cmq + Cm<x*“ To evaluate the damping in 
pitch parameter Cm^ + C m ^, semilog plots were made of the peak values 
of pitching velocity as a function of time, disregarding the algebraic 
signs of the values, and using values only from the regions where the 
longitudinal control was held fixed. From these plots, the values of 
the time required for the oscillation to damp to one-half amplitude, Tjyg, 
were obtained and applied in the equation 


In this equation) the average value of CL a over the angle-of- 
attack range covered by the particular oscillations was employed. 


Static directional stability, Cnp'.- The procedure for determining 
the values of Cnp ' was similar to that described previously for 
determining Cm^, differing only in that the sideslip-angle records were 
used to estimate the periods of the short-period oscillations. In order 
to permit direct comparison with the value of Cm a for the wing-off 
model, the coefficient Cnp* is based on the wing mean aerodynamic 
chord instead of the wing span which is used conventionally. 



Control effectiveness parameter, Cmg.- Values of C m g were deter- 


mined from the expression 


<dC m c)C m 

cd6 <dcx cdS 
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APPENDIX B 
ESTIMATION OF Cm + 

The methods used to estimate the values of the parameter C m + 
are described in the following paragraphs. 


Horizontal Tail 


The contribution of the horizontal tail to the value of C m + Cj^ 
was estimated using the expression 


Cm q + Cm^ 


2 


^Lt 


^t 


H 2 St 

"W 



+ 



The variation of with Mach number assumed for the cal- 

culations was estimated from available wind-tunnel data and is shown in 
figure 13. The value of r\^ was assumed to be 1.0 for the wing-off 

model and 0.9 for the wing-on models. The factor 1 «f was consid- 

oa 

ered equal to 1.0 for the wing-off model while, for the wing-on models, 
the variation with Mach number shown in figure 13 was obtained by inter- 
polating data from references 6 and 7* 


Fuselage 


The method of reference 8 was used to estimate the contribution of 
the fuselage to the value of C m + Cm^. Using this method, it was 
found that the fuselage contribution was reasonably close to the rough 
figure used by other investigators (reference 2) of 25 percent of the 
tail contribution. 


Wing 


There appears to be no good theory available for predicting the 

contribution of the wing alone to the value of Cm n + Cm~ at transonic 

q ex 

speeds . The trends of available supersonic theory indicate that the 
values would be quite small at the Mach numbers covered in these tests, 
and this contribution was, therefore, ignored in the estimations. 
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TABLE I.- DIMENSIONS OF FREE-FALL MODELS 1 


Gross weight, pounds 2 1238 to 1356 

Moment of inertia, slug-feet squared 2 570 to 595 

Center of gravity 2 Station 100.2 or 101.3 

Wing 

Area, square feet 9.0 

Aspect ratio 6.0 

Taper ratio 0.5 

Sweepback, quarter-chord l’ine, degrees 1+5 

Span, feet 7*33 

Mean aerodynamic chord, feet . . 1.27 

Root chord, feet 1.64 

Tip chord, feet 0.82 

Airfoil section, perpendicular to quarter-chord line 

Plane wing NACA 64A010 

Cambered and twisted wing NACA 64A810, a=0.8 

(modified) 

Twist, degrees washout at tip, measured streamwise 

Plane wing 0 

Cambered and twisted wing 10 

Incidence, degrees +0.2 

Horizontal tail (all-movable, pivoting about axis 
perpendicular to plane of symmetry of model) 

Area (including 2.0 square feet included in fuselage), 

square feet 6.0 

Aspect ratio 4.3 

Taper ratio ‘ 0.20 

Sweepback, quarter-chord line, degrees 45 

Span, feet 5*12 

1 Except as noted in the table, the wings of the two winged models were 

similar. The fuselage-tail configuration of all the models were the 
same . 

2 Two bodies of identical shape but different weight distribution were 

used. The plane-wing panels and the cambered- and twisted-wing panels 
were interchanged between the bodies during the test program so that 
the different weights and center-of-gravity positions are not readily 
identifiable with particular results. 
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TABLE I.- CONTINUED 


Horizontal tail (continued) 

Mean aerqdynamic chord (including area included 

in fuselage), feet 

Leading edge of mean aerodynamic chord Station 153.6 

Root chord, feet 1.96 

Tip chord, feet o!hO 

Airfoil section, parallel to stream NACA 65OO6 


Gap between tail and fuselage at 0 W deflection, 
inches 


Vertical tail (all-movable differentially, pivoting about 
axis perpendicular to longitudinal axis of model) 

Area (including 1.4 square feet included in fuselage), 

square feet 

Aspect ratio 

Taper ratio 

Sweepback, quarter-chord line, degrees 

Span, feet 

Mean aerodynamic chord (including area included in 

in fuselage), feet 

Leading edge of mean aerodynamic chord Station 

Root chord, feet 

Tip chord, feet 

Airfoil section, perpendicular to quarter-chord line . . . NACA 
Gap between tail and fuselage at 0° deflection. 


l/l6 


inches 


Fuselage 


Fineness ratio 

Ordinate at station x (x = 8.0 to x = 139.4), 

inches r = 8.5 


1 - 


102 


3.3 
5.0 
0.22 
• 45 
4.05 

0.93 

151.0 

1.34 

0.29 

65009 

l/l6 


12.4 

1 3/4 
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x = 8 to x=/394 Station 



Note: 

Ait dimensions are in inches. 



Figure /.- Three -view drawing of free-fall model. 
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Vertical and trans- 
verse acceleration 

g's 


Tail deflection 

deg 


Angle of attack, a , 
and sideslip angle, /9, 
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Figure 3.- Time history showing typical data obtained 
during oscillations of mode I. 
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Figure 4.- Variation of lift coefficient with angle of attack 
for the test configurations at different Mach numbers and 
horizontal- tail deflections. 
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Figure 6- Variation with Mach number of trim angles of attack 
for various horizontal - tail deflections for the test 
configurations. 
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(c) Cambered and twisted wing. 

Figure 7. - Variation with Mach number of static stability 
SC 

parameter, , for the test configurations as determined 

from fixed - control oscillations. Center of gravity at 
0.24 c. 
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Figure 8. - Variation with Mach number of aerodynamic center position of complete mode/. 
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Figure 9.- Variation with angle of attack of the pitching -moment coefficients of the 

exposed wing pane/s for the wing -on models . 
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Mach number, M 
(c) Cambered and twisted wing. 

Figure 10- Variation with Mach number of damping factor 
(C mq ■+ C m g ) for three models. 
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Figure //.- Variation with Mach number of horizontal— tail— effectiveness 
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Figure 12 % - Variation with Mach number of the static stability parameter , Cp^' t 
for the test configurations. (C n g based on SZ). 
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Figure 13,- Characteristics of models assumed in 
estimating values of (C m ^ + C m ^ ). 


NACA-Langley - 8-29-52 - 400 


CONFIDENTIAL 




SECURITY INFORMATION 


CONFIDENTIAL 





CONFIDENTIAL 


